On the density-potential mapping in time-dependent density-functional theory 
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Abstract 

The key questions of uniqueness and existence in time-dependent density-functional theory are usually formulated only for 
potentials and densities that are analytic in time. Simple examples, standard in quantum mechanics, lead however to non- 
analyticities. We reformulate these questions in terms of a non-linear Schrodinger equation with a potential that depends 
non-locally on the wavefunction. 
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I. INTRODUCTION 

The reduced one-particle probability density (hence- 
forth density tout court) in a gas of identical interact- 
ing particles in an external confining potential carries a 
wealth of information. For example, the linear response 
of the density to a driving field is related to density- 
density correlations and the excitation spectrum of the 
system. By continuity, the density evolution sets the to- 
tal current through an arbitrary open surface across a 
finite system driven out of equilibrium, giving access to 
transport properties, or across a closed surface, giving 
access to transfer processes out of a chosen volume. The 
density determines the generalized mean forces on macro- 
scopic degrees of freedom, with which the gas interacts, 
and hence the coupled dynamics of the two subsystems. 

In fact, the theorems of time-dependent (TD) density- 
functional theory (DFT) prove that, given an initial 
many-body state, the TD density contains in principle 
the expectations of all observables of the interacting sys- 
tem, evolving under a TD scalar potential. Following 
(the much older) ground-state DFT [l|,[|, TDDFT calcu- 
lations operate by constructing a fictitious system of non- 
interacting electrons, evolving under another TD scalar 
potential, such that the fictitious system reproduces the 
electron density of the interacting system. TDDFT thus 
provides access to the properties of the interacting sys- 
tem via the much more tractable corresponding non- 
interacting problem. In recent years, TDDFT has be- 
come a central tool for a range of problems involving 
departures from the electronic ground state (GS). Most 
applications are in the linear-response regime, calculating 



the electronic excitation spectra and response of atoms, 
molecules, solids, even biomolecules, exposed to external 
fields Real-time electron dynamics in strong fields 
have also been studied, including electron transfer pro- 
cesses and transport, ionization, high-harmonic genera- 
tion^and coupled non-adiabatic electron-nuclear dynam- 
ics 



TDDFT is based on two distinct results. The first 
is that, under a given particle interaction W and for a 
given initial many-body state |^(0)), there is a 1:1 corre- 
spondence between the ensuing evolution of the density 
and the external potential acting on the system. This 
is the Runge-Gross theorem (RG) Q. It implies that 
the potential and all other properties of the system are 
functionals of the density and initial state, and practical 
calculations hope to approximate some of these function- 
als accurately. 

The second result is that for every density n(r,t), 
evolving from a given initial state 1^(0)) under a given 
particle interaction W and external potential v(r,t), 
there exists an external potential v'(r, t) that returns the 
same density, n(r,t), under another particle interaction 
W' (the key case of interest being W' — 0), starting 
from, in general, another initial state ^'(O)) (so long as 
the two initial states share the same n(r, 0) and n(r, 0)). 
This result is van Leeuwen's u-representability theorem 
(vL) [6j , later extended to TD current DFT 7} . It means 
that, by contrast with GS DFT, there is no doubt about 
the existence of a KS system, i.e. a set of non-interacting 
electrons whose TD density matches that of the interact- 
ing system. The w-representability difficulty is lessened 
in the TD case because the TD Schrodinger equation 
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(TDSE) is first-order in time, and the dependence on the 
initial wavefunction takes care of a large part of the dif- 
ficulties associated with v-representability. 

The first result has been proven for potentials that are 
analytic in time about the initial time (henceforth de- 
noted t-analytic), i.e. such that v(r, t) is equal to its 
Taylor series expansion in t about the initial time, for a 
finite time interval. The second assumes i-analytic po- 
tentials and densities. There are also two extensions of 
RG in the linear-response regime thatgo beyond these 
analyticity requirements. In the first [8j, the short-time 
density response to "small" but arbitrary potentials has 
been shown to be unique under two assumptions: that 
the system starts from a stationary state (not necessar- 
ily the GS) of the initial Hamiltonian and that the cor- 
responding linear density-response function is t-analytic. 
In the second [9], uniqueness of the linear density re- 
sponse, starting from the electronic GS, was proven for 
any Laplace-transformable (in time) potential. As most 
physical potentials have finite Laplace transforms, this 
represents a significant widening of the class of poten- 
tials for which a 1:1 mapping can be established in the 
linear- response regime, from an initial GS; this includes 
e.g. potentials turning on as e~ c / tn with C > 0, n > 0, 
or t p with p positive non-integer. It covers most cases of 
physical interest, under these conditions. 

The TDSE is a coupled partial differential equation, 
which inextricably links analyticity properties in space 
and time. Thus potentials and initial wave functions that 
are (or could be) standard examples in textbooks, in 
which the potential or the initial wavefunction has some 
non-analytic properties in space, yield densities that are 
not analytic in time. We demonstrate this with several 
examples in the last part of this paper. This implies that 
the assumptions underlying vL are more restrictive than 
they may first appear. 

In this paper we build on the foundational work in 
Refs. p5, 6] on the problems of uniqueness and existence 
in TDDFT. Our goal is to extend the fundamental frame- 
work laid there, by suggesting a simultaneous formulation 
of these two problems, which overcomes the requirement 
that the potential or density be analytic in time. This 
formulation is close to (but different from) several prob- 
lems whose solutions have been proven to exist mathe- 
matically. We relate the one-to-one mappings to the ex- 
istence and uniqueness of solutions to a particular type of 
non- linear TDSE (NLSE) . While there has been a consid- 
erable amount of work on the NLSE licj, the particular 
non-linear structure we have here has not been investi- 
gated before, to the best of our knowledge. We do not 
give a general solution to this problem, but discuss the 
features of the equation in relation to other NLSEs in the 
literature. As the existence of solutions of the latter do 
not require ^-analyticity, this suggests that the reformu- 
lation will be a useful tool for exploring the limitations 
(or extensions) of our knowledge of the domain of the 
Runge-Gross functional. 

In Section II, we carefully introduce our reformulation. 



We explore some consequences of the reformulation and 
what is known about solutions to the NLSE that we find. 
Section III turns to the motivation behind this endeavor, 
with the aid of simple examples where the density evolu- 
tion under t-analytic potentials is non-analytic in time. 

II. ALTERNATIVE FORMULATION OF THE 
DENSITY-POTENTIAL MAPPING 

It is straightforward to see that, for a given initial 
state 1^(0)) for a system of N electrons, a given state- 
evolution, \^(t)), can be generated by at most one TD 
scalar potential v(r,t), as follows from TDSE 

m i*(t)> = H(t) i*(t)> = (f+w+ op)) m)) (i) 

as long as the wavefunction ^(ri, rjv, t) is non- 
vanishing on a dense set of values of its arguments. Here, 
T, W and v(t) are the many-electron operators for the to- 
tal kinetic energy, the electron-electron interaction (with 
W = corresponding to non-interacting electrons) and 
v(r, t), respectively. That no two distinct #'s can generate 
the same up to a physically irrelevant purely time- 

dependent phase, can be shown by assuming otherwise, 
inserting into TDSE, and finding t)i(t) — i>2(t) — c(t), a 
physically irrelevant purely time-dependent constant. 

What RG proves is the much stronger statement that 
not only the states must be different, but also their one- 
particle densities, 

n(r,t) = (*(t)|fi(r)|tf(t))=n*(r,i), (2) 

must differ, under the restriction that the potentials 
considered are analytic around t = 0. Here n(r) = 

Y^iLi l K r — *i) i s the number-density operator for the 
system of N electrons with position operators {f^}. This 
restriction means that potentials that switch on at t > 
like e _c /* , t p with p non-integer, or t n ln(i) are not cov- 
ered by the proof. (Note that the first example is in- 
finitely differentiable, with vanishing derivatives at t = 0, 
while higher-order derivatives of the second and third 
types diverge as t — > + ). These potentials are Laplace- 
transformable in t, and in the linear-response regime 
starting from the GS are covered by the vL extension 
of RG Q. 

Here, we formulate (but do not prove) a 1:1 mapping 
defined via the existence and uniqueness properties of 
the solution to a system of coupled equations. Differen- 
tiating the continuity equation (n(r, t) + V • j(r, t) = 0, 
where j(r, t) is the current density) once with respect to 
t gives [6| 

— V-[n*(r,i)Vu(r,t)] = n*(r,i)+V-a[*(t)](r,t) , (3) 
m 

where 

a[tt(t)](r,t) = ± (*(t)|[j(r),T + W]\9(t)) (4) 
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with j(r) the current-density operator 

We now choose a given prescribed density, which we 
denote by n a i m (r,i), and a given initial many-electron 
state |^(0)), such that n*(r, 0) = n a ; m (r,0), ?W(r, 0) = 
n a i m (r, 0). We then ask, does a potential v(r,i) exist, 
such that the density it generates via Eqs. (p} and 
([2]), from the chosen initial state, is equal to the pre- 
scribed density, u*(r, t) — n a j m (r, t), and is this potential 
unique? 

Noting the appearance of the time-evolving state |^(t)) 
in the evaluation of a[$(t)](r,t) in Eq. @, we re-phrase 
the question by asking whether the coupled equations 



1 (r,i)+V-a[*(i)](r,t) 



ift|*(t)) = H(t)\*{t)) 



(7) 
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V • [n aim (r,t)Vv(r,t)] = n aim (r, t) + V • a[tf(i)](r(3) 



ift|*(t))=ff(t)|*(t)) 



(6) 



have a unique solution for u(r, t) and |vP(f)), given the 
inputs ^(O)) and n a j m (r, t). At a notional level, these 
coupled equations work as follows. Eq. (0 determines 
the potential i>(r, t), at every t, as an instantaneous func- 
tional of the state the (fixed-time) equation is of 
Sturm-Liouville form, for which a unique locally square- 
integrable solution exists for strictly positive densities, 
and locally square-integrable right-hand side [l2j. The 
restriction to nodeless densities is in practice a minor 
one as most densities of iV-electron systems do not have 
nodes; but it certainly excludes a small subset of them, 
and excludes typical TD and excited-state densities of 
one-electron systems. The restriction to locally square- 
integrable right-hand sides (meaning that the integral of 
the square of the right-hand side over a finite region is fi- 
nite) is generally satisfied for physical wavefunctions [l3[ . 
A formal solution to Eq. (|5|) is discussed in the appendix. 
Eq. © then determines the evolution of the state, un- 
der the potential just obtained, through an infinitesimal 
time-step to the "next" time, and so on. Eqs. ([5]) and 
© therefore define the density-potential mapping: if a 
unique solution (w(r, t), \^(t)}) to Eqs. ([5]) and ^ exists 
for a given (|^(0)), n a i m (r, t)), then there is a 1:1 map- 
ping between the density and potential. If, further, the 
density n a i m (r,i) is w-representable, it directly follows 
that the unique solution (v(r,t),\ty(t))) reproduces this 
target density, n#(r, t) = n a i m (r,t). 

Non-u-representability can arise in two ways. First, 
if no solution exists to Eqs. ([S])-® for a particular 
"aim(r, t), then ri a i m (r,t) is not v-representable. Second, 
a solution to systems of equations of type ©-(El) does not 
necessarily have n^(r, t) ^ n aim (r, t) jl4| . If there does 
exist a unique solution (v(r, t), \^(t))) but one that fails 
to return the chosen target density, n*(r, t) ^ n a ; m (r, t), 
then this target density is not w-representable [la |. 

The coupled equations ([5]) and ^ can be expressed as 
evolution equations as follows. By taking one further t- 
derivative of Eq. ([5]), we obtain the system of first-order, 
in time, partial differential equations 

— V- [n aim (r,i)Vu(r,t)] = — V • [n aim (r, t)Vv(r, *)] 
to m 



for the pair (v(r, t), \ty(t))). The initial conditions are the 
initial state |W(0)) and v(v, 0), given by the solution of 

- V • [n aim (r, 0)Vu(r, 0)] = n aim (r, 0) + V • a[*(0)](r, 0) . 

TO 

(8) 

A system of first-order ordinary differential equations 
with fixed initial conditions always has a unique solution 
(given appropriate continuity conditions on the right fl6l | , 
such as Lipschitz continuity). But partial differential 
equations are more complex. While the initial-value 
problem for many classes of partial differential equations 
is known to have a unique solution, to our knowledge, 
there are no results yet for equations of the form ([7)1, 

Instead, we consider a reformulation of the coupled 
equations ([5]) and © as a NLSE with non-local potential 
term, i.e. the potential at r depends on the wavefunction 
not just at or near r but further away too. The solu- 
tion for v(r, t) is an instantaneous functional of ri a j m (r, t), 
?iaim(r, f) and |\&(t)) as can be seen from Eq. (T5J), and its 
dependence on \^(t)) is quadratic. Denoting this solution 
symbolically by v{r.t) = v{n aim , n aim , *](r, t), Eqs. © 
and ([6]) now reduce to the effective NLSE 

ffi|*(t)) = (f + W + v[h aim , naim ,y](t)) |*(t)). (9) 

The problems of u-representability of the chosen density 
7"iaim( r )i) and of the uniqueness of this representability 
reduce to the problem of the existence and uniqueness of 
the solution to Eq. © (together with the requirement 
that the solution returns nq,(r,t) — n aim (r,t), whose 
negation implies that the chosen density n a i m (r,£) is not 
w-representable, as discussed earlier). Two specific issues 
that bear on this problem are the nature of the non- 
linearity of this equation and the boundary behaviour of 

*](*)■ 

NLSEs with non-local interaction were investigated by 
Ginibre and Velo [l7j . who in particular proved under 
certain conditions the existence and uniqueness of short- 
time solutions. However their non-linear term (a) does 
not include Coulomb interactions, (b) deals only with 
the version of Eq. (|5|) for the ordinary Laplacian, and (c) 
does not involve derivatives of the wavefunction. Much 
of the recent mathematical work involves non-linearities 
depending on the wavefunction and a single derivative, 
whereas in this case the potential depends non-locally 
on four derivatives of For point (b), the results of 
Ref. 17[ could be extended to a general uniformly el- 
liptic operator on the left (i.e. one which has eigenval- 
ues bounded below by a positive constant, everywhere in 
space), but this is not what we have: the density appear- 
ing on the left-hand side of Eq. (|5|) decays at infinity, so 
the differential operator on the left is not uniformly el- 
liptic (i.e. n a i m (r, t) is not gapped away from zero). This 
makes it difficult to obtain estimates on the behavior of 
w(r, t) at infinity needed to reproduce the Picard iteration 



3 



needed to show existence and uniqueness as in Ref. (171 ]. 
Problem (b) is absent in the simplified case of periodic 
boundary conditions, and smoothing over Coulomb sin- 
gularities deals with problem (a) (as would be relevant 
for solid-state calculations using pseudo-potentials with 
a smoothed-out interaction term). For example, Burq, 
Gerard, and Tzvetkov [HI prove short-time existence and 
uniqueness of the NLSE with local potential term on com- 
pact manifolds. However, even in this case it is not clear 
whether the given NLSE has the hoped-for solution, be- 
cause of the loss of derivatives (point (c) above) . Ref. [l9[ 
proves well-posedness for a NLSE that loses two or more 
derivatives, however for Hamiltonians with local interac- 
tions only. Part of the motivation for this paper is the 
hope that the further analysis of this problem in mathe- 
matical physics will stimulate future work. 



III. EXAMPLES OF NON T-ANALYTIC DENSI- 
TIES 



The practical motivation behind a formulation that 
goes beyond the analyticity requirements of the origi- 
nal Runge- Gross and van Leeuwen theorems is that it is 
not unusual, and certainly not pathological, for temporal 
non-analyticity to arise. If it were true that under any 
particle interaction, a i-analytic potential always gener- 
ates a t-analytic density, and vice-versa, then the con- 
structs of Refs. [5} and p] would be sufficient to provide 
a self-contained (even if somewhat restrictive) framework 
for TDDFT. However, we now argue, with the aid of sev- 
eral examples, that in general a i-analytic potential need 
not produce a i-analytic density. Although the RG 1:1 
mapping holds in such cases, the u-representability proof 
of Ref. [6| does not apply, opening up the question of the 
existence and uniqueness of a corresponding Kohn-Sham 
system. 



A. Spatially non-analytic initial states 

In the first set of examples, spatial non-analyticity in 
the initial wavefunction leads to temporal non-analyticity 
in the density. 



1. "Bump" initial states: smooth but non-analytic 

Consider a free particle in Id with an initial wavefunc- 
tion tj)(x, 0) that is constant in space within a region, e.g. 



f(x) for x < a 
4>{x, 0) = { C for a < x < b 
g(x) for x > b . 



(10) 



Here, C is a constant, and f(x) and g(x) are assumed 
to be such that ip(x, 0) is infinitely differentiable in x 



everywhere. By repeated application of TDSE 



in^{x,t) = -—i/j"{x,t), 
Zm 



(11) 



we see that, for a point x = Xq in (a, b), all i-derivatives 
of ip(xo,t) vanish at t = [20] ■ Thus, the Taylor series 
in time for ij>(xo, t) at fixed xq is 



Mxo,t) = J2 ^dfy(x ,t)\ t=0 = C. 



fe=0 



(12) 



At the points xq considered, this series exists, and con- 
verges (with, furthermore, an infinite radius of conver- 
gence), but not to the function itself: ipr(xo,t) = C ^ 
^(xo, t). The true wavefunction, ip(xo, i), evolves in time, 
since the initial state is not a free-particle eigenstate. In 
particular, the true solution will have density from the 
outer regions, x < a and x > b, dispersing into the inner 
region, a < x$ < b. Therefore, ip(xo,t) is not analytic in 
time. Similarly, using 



d?n(x ,t)\ t=0 = J2 



{k-l)W. 



x dt l r-(xo,t)\t=od l Mx ,t)\ t=0 (13) 

for the Taylor series for the density at xq we obtain 
n,T(xo,t) — |C| 2 ; it exists and converges, but not to 
n{xQ,t): riT(xo,t) ^ n(xo,t). Thus, the density is not 
i-analytic, even though the potential (v(x, t) = 0) is. 

The above arguments may be straightforwardly ex- 
tended to initial states that have compact support but 
are smooth: there in the region (a, 6), ip(x, 0) = f(x) 
while outside this region, ip(x, 0) = 0. Such functions are 
sometimes referred to as "bump functions" . The true so- 
lution i/j(x,t) disperses, while the time Taylor series of 
ip(x,t) does not. Although formally solving the TDSE, 
the Taylor series is not a valid solution because it is not 
uniformly convergent in the non-constant region. Simi- 
lar observations and a detailed analysis for free propaga- 
tion of smooth compactly supported initial states can be 
found in Ref. [HI. 

The above example is easily extended to any initial 
"spline" state 



■ f{x) for x < a 
ip(x,0) = < i} n (x) for a < x < b 
g(x) for x > b , 



(14) 



where ipn(x) is the nth stationary state, with energy E n , 
of some static confining potential v(x), and f(x) and g(x) 
are arbitrary (but may, if desired, be assumed to be such 
as to ensure infinite smoothness everywhere). Starting 
from this initial state, we consider evolution under the 
potential v(x). By repeated application of TDSE 



ihip(x, t) 



2m 



ip"{x,t) +v(x)tp(x,t), 



(15) 
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for Xo € (a, b) we arrive at iPt(xo, t) — e i/j n (xo) = 

e- i£ Wft^( XO) o) and nr(x ,t) = IV^o)! 2 = n(a;o,0), 
with the conclusion, again, that neither the wavefunction 
nor the density is i-analytic. 



2. Initial states with a cusp 

Consider evolving an initial state with a cusp in free 
space. We can solve this explicitly in Id, with: 



ip(x,0) = ^e"^ 1 
Using the free-particle propagator 



and the integral representation 



^i(x-x') m/2ht 



ip{x,t) 
we obtain 

il>(x,t) = ^ 



G(x, t; x , 0)ip(x , 0) dx 



(16) 



(17) 



(18) 



.^\tiK t/2m J pKi 



1 — erf 



x + ihnt/m 
y/2ifit/m 



1-erf 



-x + ihnt/m 
m 



(19) 



Immediately from t = + , the cusp vanishes, and the ini- 
tial exponential gradually disperses and develops oscilla- 
tions in the tail (see Fig. [TJ and note the diminishing 
scale on the vertical axis). 



clearly not equal to the true wavefunction ip(x,t) 
(Eq. {T5J}). Noting that the true solution to TDSE is 
defined and continuous everywhere, we define ipT{0,t) 
such as to make ipT(x,t) continuous. But then it obeys 
the TDSE 

ihip T (x,t) = -^-^(x,t) - —5(x)i> T (x,t) , (21) 
Ira m 

as opposed to the free-particle TDSE, Eq. (fTI) . obeyed by 
ip(x, t). Therefore the Taylor series gives a delta- function 
error in solving the TDSE. 

The electron-nuclear Coulomb singularity in 3d gives 
rise to cusps in stationary states and densities in real sys- 
tems, i.e. initial states are typically not smooth at the 
singularity. For example, the GS of the hydrogen atom, 

iPhM = w^-re -r / a °, where ao is the Bohr radius, has 

V 7ra 

analogous mathematical structure to the initial state of 
Eq. (1161) . If made to evolve in a potential in which the 
Coulomb singularity changes (e.g. a nucleus that van- 
ishes, admittedly unphysical), the space-time coupling in 
the TDSE results in a non-analytic evolving density, in 
a way that closely follows the above Id model. We need 
not however think of the initial state as a GS of any po- 
tential: the initial state in Eq. (ffjj)) or its 3d analog iI>h(t) 
(or any of the examples given here) can be prepared in 
many different ways, for example, using particular inter- 
ferences of free-particle eigenstates. The essential point 
is that we have a quantum-mechanically allowed initial 
state, which, when freely-evolved, gives a non-t-analytic 
density. 

That non-t-analytic behavior is not a consequence of 
a change or removal of a singularity in the potential, is 
shown explicitly in the next example. 
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FIG. 1: The density of an initial exponential dispersing under 
free-particle evolution. 

From the Taylor series, away from x — 0, we obtain 
simply 



i/j t (x ^ 0, t) = v^e mK2 */ 2m e- KN 



(20) 



3. Free evolution of the GS of a smooth bounded potential 

Consider a free particle in Id evolving from 

V>(x,0) =M{l-e~ c ' x2 ), (22) 
where C > is a (real) constant and J\f = 



1/ Y 2-\/27rC(v / 2 — 1) is the normalisation constant. This 
initial wavefunction can be thought of as an example of 
the earlier bump states, with the interval [a, b] shrunk to a 

point. It is the GS of the potential v(x) — — , i C \ , 

infinitely diffcrentiable; see Fig. 2 for a plot. 
Using Eqs. JT7J and (JTHJ), we find 



ip(o,t) =j\f (i - e -V 2mC / mt ) 



(23) 



and 



»(0, t)=N 2 (l - 2 e-V m C/ftt CQS ^ mC / ht + e -2y/mC/ht^ 

(24) ' 
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FIG. 2: The initial state of example A. 3 (taking C = 1), and 
the potential in which it is the ground state. 

Neither is t-analytic. The non-analyticity in the density 
takes the form of temporal evolution that picks up in- 
finitely slowly. 

To summarize, the examples above illustrate how a 
non-t-analytic density may arise from a (trivially) t- 
analytic potential, starting from an initial state that is 
spatially non-analytic, even though it could be smooth. 
Indeed Ref. points out that even for the simplest 
free-particle evolution, starting from an initially smooth 
state (our examples 1 and 3 are examples thereof), al- 
though the Taylor series solution is a formal solution of 
the TDSE, it is a valid solution only for a very restricted 
class of initial states, due to non- convergence. Even the 
time-evolution of the density in a hydrogen atom in an 
electric field can be shown to be non-i-analytic [22| . The 
Runge- Gross theorem applies to all cases, as it is a refor- 
mulation of the quantum dynamics of any initial state 
allowed in quantum mechanics (satisfying appropriate 
decay at infinity) under a t-analytic potential. But in 
cases where the ensuing evolution is non-i-analytic, the 
w-representability theorem does not apply. 



B. Spatially non-analytic potentials 

Finally, we consider potentials with the spatial non- 
analyticity of "bump functions" . We have one particle in 
Id, starting from the nth stationary state i/j n (x), with en- 
ergy E n , of some smooth static confining potential vq{x). 
Now evolve ift(x, t) (from ip{ x , 0) = ip n (x)) under the po- 
tential 

v(x) = V (x) + Uextfx) , (25) 

where 

!f(x) for x < a 
C for a < x < b (26) 
g(x) for x > b . 

Here, C is a constant, and f(x) and g(x) may be chosen 
such as to give infinite differentiability in x everywhere. 
By applying TDSE repeatedly, for xq £ (a, b) we arrive at 



and riT{xo,t) = \ip n {xo)\ 2 — n(xo,0). Since both the 
wavefunction and the density must evolve, we conclude 
again that neither of them is ^-analytic, even though the 
potential is. 



IV. SUMMARY AND OUTLOOK 

TDDFT is an exact reformulation of quantum mechan- 
ics in which the density replaces the wavefunction as the 
basic variable. The density-potential mapping is proven 
for any initial state evolving under a i-analytic potential, 
such that the density decays appropriately at infinity. 
The mapping to the KS system is proven for t-analytic 
densities and potentials. However, quite generally, the 
densities arising from the TDSE are not ^-analytic, as 
the simple examples in Sec. IIIII demonstrate. The initial 
states and potentials considered are certainly within the 
realm of traditional quantum mechanics, and nothing in 
the Runge- Gross theorem rules them out from consider- 
ation for TDDFT. 

The examples suggest that in general a ^-analytic po- 
tential does not give a i-analytic density. However, then 
it is no longer clear that a density, generated by a t- 
analytic potential under one particle interaction, is rep- 
resented by a i-analytic potential under another. This 
would throw into question the existence and uniqueness 
of the non-interacting KS system. Therefore, an exten- 
sion of the TDDFT framework, beyond i-analytic poten- 
tials and densities, is necessary. The present proposal is 
that Eqs. (J5J) and ([5]) (or their two further formulations 
in Section |H| provide a route to one such extension. 

The framing of the fundamental theorems of time- 
dependent density-functional theories in terms of well- 
posedness of a type of NLSE first appeared in Ref. [23[ 
where it arises naturally in Tokatly's Lagrangian formula- 
tion of TD current-DFT, known as TD-deformation func- 
tional theory. The traditional density-potential mapping 
question is avoided in TD-deformation functional theory, 
where instead this issue is hidden in the existence and 
uniqueness of a NLSE involving the metric tensor defin- 
ing the co-moving frame. Investigating the relation be- 
tween the particular NLSE of Ref. [23| and the one in the 
present paper is an interesting avenue for future work. 

It is hoped that these brief considerations will moti- 
vate lively work to take these preliminary developments 
further. 
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VI. APPENDIX 

In this appendix, we explore some interesting proper- 
ties of the solution to Eq. ((5]). Specifically, we discuss 
a reformulation that makes contact with the more stan- 
dard theory of second-order elliptic operators, namely, 
operators that are uniformly elliptic at infinity. 

We can formally solve Eq. ([5]) for v(r,t), with the aid 
of the transformation 



w(r,t) =v(T,t) y /n aim (r,t). (27) 
Then Eq. §5§ becomes 

A(t)w(r, t) = -to / S(r,t) = = u(r, t) , (28) 



where 



y/n aim (r,t) 



1 V 2 n aim (r,t) 1 , 

* "'aim M) 



|29) 



A ^ V + 2 n aim (r,t) 4 V n aim (r,t) 
«(r,t)=n ailn (r,t) + V-a[*(t)](r,t). (30) 

We must now solve this equation for ti)(r, t). In the fol- 
lowing, t may be treated as a parameter. The operator 
A(i) has the structure of an ordinary one-particle Hamil- 
tonian. <po(r, t) = a/ n a i m (r, t) is an eigenfunction of A(t) 
with eigenvalue 0. Since f n a j m (r, t)dr — N, <f>o(r,t) is 
normaliseable (and thus obeys the boundary condition 
4>o(r,t) — > as |r| — > oo). Therefore, cj>o(r,t) is a bound 
state of A(t). Since the eigenvalue of </>o(r, i) is zero, for 
a solution to Eq. (|2"8"]l for w(r,£) to exist, it is necessary 
that the right-hand side be orthogonal to (f>o, 



<Mr,i) 



*(r,t) 



\/naimO 



= dr = J s(r,t)dx = 0. (31) 



Via Eq. ([30]). this requires J g a[*(i)](r,i) • dS to vanish 
as the surface S expands to infinity, which is a condition 



on the state in the solution of Eqs. © (or ([28])) 

and ([6]). Since the eigenvalue of <fo(r, t) is zero, further- 
more, Eq. ([28} defines tu(r, i) only up to an arbitrary 
additive amount of <fio(r,t). We choose this amount to 
be zero. Denoting the eigenfunctions and eigenvalues of 
A(t) by {<j) v {T,t)} and {A„(t)}, respectively, we can then 
uniquely solve Eq. ([28]) for w(r,t): 



w(r,t) 



J g{v,v',t)u{v',t)dv' , 



where 



(32) 



(33) 



The potential v(r,t) is then obtained from Eq. ([2"T]) . The 
given solution for w(r, t) by construction is orthogonal 
to cf>o(r,t), and corresponds to a specific choice of the 
additive spatial constant in v(r,t) [24|. The boundary 
behaviour at infinity of w(r,t) in Eq. ([32[) . and hence 
of v(r,t) in Eq. (j27|). is determined by Q{t) and by the 
quantity u(r,t) in Eq. (|2"8"]). 

No difficulty with this solution for w(r, t) arises, if the 
target density n a ; m (r, t) is nodeless. Then <fi (r,t) too is 
nodeless, and since it is a bound state of A(t), it must be 
the GS. Then all other eigenvalues of A(t) are positive, 
Q{t) in Eq. (|3"51 exists, and the potential u(r,t) in Eq. 
(f27|) is defined everywhere. 

The situation is more complicated if n a ; m (r,i) has 
nodes. Then, first, A(t) may have bound states degener- 
ate with 0o ( r j t) an d the requirement of Eq. pip must be 
extended to all such states, which constitutes a stronger 
condition on s(r, t). (That granted, the formal solution 
for w(r,t) is then still given by Eq. (|3"2"j) . but with all 
bound states of A(t), degenerate with (po(r, t), excluded 
from the definition of Q{t) in Eq. (13"3"1) .) Second, both 
u(r,t) in Eq. ([28]) and v(r, t) in Eq. (1271) may diverge 
at the nodes of n a im(r, i). One would then have to con- 
sider the nature of the singularity, and its effect on the 
propagating state vector. 
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